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£N) , The field-angle dependent specific heat and thermal conductivity in the vortex phase of UPd2 AI3 is 

studied using the Doppler shift approximation for the low energy quasiparticle excitations. We first 
give a concise presentation of the calculation procedure of magnetothermal properties with vortex 
and FS averages performed numerically. The comparison of calculated field-angle oscillations and 
the experimental results obtained previously leads to a strong reduction of the possible SC candidate 
states in UPd2Al3. The possible SC gap functions have node lines in hexagonal symmetry planes 
containing either the zone center or the AF zone boundary along c. Node lines in non-symmetry 
planes can be excluded. We also calculate the field and temperature dependence of field-angular 
oscillation amplitudes. We show that the observed nonmonotonic field dependence and sign reversal 
O ' of the oscillation amplitude is due to small deviations from unitary scattering. 
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I. INTRODUCTION 



The U-based heavy fermion (HF) superconductors (SC) are supposed to have unconventional SC order parameters 
usually (but not necessarily) associated with anisotropic gap functions A(k) that have node points or lines on the 
Fermi surface (FS) where A(k) = [l|,|2||. This is thought to be the result of a purely electronic mechanism of Cooper 
pair formation which favors anisotropic gap functions due to a strong on-site heavy quasiparticle repulsion. The 
type and position of nodes in A(k) is intimately connected with the symmetry class to which A(k) belongs, in most 
cases described by a single irreducible representation of the high temperature symmetry group 0, Q . Experimental 
evidence for the presence of node lines is obtained from thermodynamic and transport quantities as well as resonance 
experiments, but to locate their exact position on the FS and hence restrict the number of possible representations 
for the SC order parameter is a difficult task. For example in UPt3 it took a considerable time until it was identified 
as the two-dimensional odd parity (spin-triplet) E2 U representation and still there is no unanimous agreement on this 
\ symmetry |5j. 

Recently the experimental determination of SC gap symmetries has been much facilitated by the advent of a new 
I/"") ! method, namely the investigation of field-angle dependence of specific heat and thermal conductivity at tempera- 
tures T <C T c . From the typical angular oscillations observed in these quantities under favorvable conditions (small 
quasiparticle scattering) one may deduce the position of the nodal lines or points of A(k) with respect to the crystal 
axis. Knowledge of these positions narrows down the possible choices of representations for A(k) considerably. This 
method has been sucessfully applied to unconventional organic SC , to ruthenates @ , borocarbides [H Q and 
Ce,Pr-based HF superconductors 0,0,^2- ^ ^ s based on the 'Volovik effect' 0] which means the appearance of 
quasiparticle states in the inter-vortex region of unconventional SC due to the presence of nodes with A(k) = along 
certain directions in k-space. The zero energy density of states (ZEDOS) of these continuum states and hence their 
O ■ contribution to specific heat and thermal conductivity depends on the relative orientation of field direction, nodal 
positions and crystal axes through the superfluid Doppler shift (DS) effect of quasiparticle energies. This results in 
• i-H . the angular oscillations of specific heat and thermal conductivity which contain important information on the nodes 
' of the gap function. 

In addition this method has now been applied for the first time to a U-based superconductor |T^. namely the 
intermetallic moderate HF (7 = 140 m ";j^ ) compound UPd2Al 3 ^5|. This compound was in the focus of interest 
in recent years because it is the only HF superconductor where direct evidence for the microscopic nature of the SC 
pairing mechanism has been found. This is connected with the fact that UPd2Al3 is the most clear cut example of 
a U-based superconductor (SC) with dual-nature 5f electrons, some of which are localised and some itinerant. The 
former can be considered as 5? CEF states and the latter as conduction band states The mass enhancement of 
conduction electrons (m*/m(, ~ 10, m& is the band mass) is a result of their interaction with the propagating CEF- 
excitations ('magnetic excitons') associated with the localised 5f electrons. The induced-moment AF order with Tn = 
14.3 K, Q = (0,0,0.5) (r.l.u.) and moderately large /z = 0.85^b coexists with SC below T c = 1.8 K. In complementary 
INS [TtI Il8| and quasiparticle tunneling experiments |19| both 5f components were investigated and it was concluded 
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18] that magnetic excitons mediate Cooper pairing. Theoretically this new pairing mechanism was investigated in 
|2Ct |2*H and possible symmetries of the SC states were discussed, also in connection with existing Knight shift [2^ 
and upper critical field results [23] ■ The conventional itinerant spin fluctuation mechanism has been investigated in 
[2II2 H and also in po | . 

The plausible SC gap functions obtained in |2l| from a microscopic model predict a node line parallel to the 
hexagonal ab-plane but several T)qh representations with different parity are possible solutions. Furthermore the 
alternative spin fluctuation model of |24j predicts node lines perpendicular to the basal plane. Therefore further 
investigation of the gap structure of UPd2Al3 has turned out to be necessary. It was already suggested in [26| that 
field-angle resolved experiments might be helpful to clarify the situation. They have now indeed been performed in 

The purpose of this paper is twofold: Firstly, although the theory of magnetothermal properties in superconductors 
on the basis of the DS approximation is well developed, the results are scattered through the literature and therefore 
we first give a concise and complete outline of the necessary computational steps for SC with uniaxial symmetry in 
the superclean limit. The calculation of linear specific heat coefficient y(T,H) = C(T,H)/T and thermal conductivity 
Ku (T,H) (i = x,z,y) involves averaging over both quasiparticle momenta and energies and the vortex coordinate. 
For quantitative predictions the five-fold integrations are carried out fully numerically for each of the candidate gap 
functions. Also this has the advantage that one can study the temperature dependence of oscillation amplitudes and 
real FS geometry effects. Secondly we want to apply the DS theory in detail to UPd2Al3 and study the predicted 
field-angle variations of the above quantities for the possible gap functions with special emphasis on the problem of 
node- line position along c* . We also discuss the influence of FS cylinder corrugation on field angle dependence and 
the temperature variation of angular oscillation amplitudes and investigate the dependence on the scattering phase 
shift. 

In Sect, [n] we introduce the concept of the Doppler shift approximation for quasiparticle energies and in Sect. ITTT1 
we give the explicit expression for this quantity in two FS geometries. In Sect. llVl we define the necessary averages 
over vortex coordinates (superfluid velocity field) in the single-vortex approximation. The calculation procedure for 
specific heat and thermal conductivity in the superclean limit is given in Sect. E| Then in Sect. E3 we apply the 
theory to UPd2Al3 and discuss the results for the most prominent candidate gap function in view of the available 
experimental results in Finally Sect. IVIII presents our conclusion on the gap symmetry of UPd2Al3 and an 

outlook on theoretical developments. 



II. DOPPLER SHIFT OF SC QUASIPARTICLES IN THE VORTEX PHASE 

In the vortex state the superfluid has acquired a velocity generated by the gradient of the condensed phase as given 
by 

2e 

v s = ftfn _1 (V(/) — — A) = ftm _1 V</> (1) 
nc 

It is connected to the screening current circulating the vortex by 

j s (r) = 2e« s (r)v s (r) = fvxh(r) (2) 

47T 

Here = m a 5ij + (m c — m a )hihj is the uniaxial mass tensor with n giving the symmetry axis. Furthermore n s (r) 
is the superfluid density and h(r) is the local magnetic field strength. The above equations hold for large Ginzburg 
parameters when local London electrodynamics is applicable |27l |28( . 

Quasiparticle excitations out of the condensate with momentum have an energy Ei(p^) in the local rest frame 
of the superfluid. The transformation to the laboratory frame is given by the universal law |29j 

P = PL 

E(p,r) = £ L (p)+p-v s (r) (3) 

The second equation may be interpreted as a Doppler shift (DS) of quasiparticle energies due to the moving condensate. 
In an unconventional superconductor with gap nodes A(k) = on the FS low energy (E <C A) quasiparticles can 
tunnel to the intervortex region where they acquire a DS according to Eq. @. Since in a nodal SC the zero-field 
DOS starts like a power law N(E) ~ E" the DS, after averging over position r and momentum p, will lead to a 
non-vanishing ZEDOS N(E=0,H) of quasiparticles which determine the low temperature (T -C T c ) specific heat 
and thermal transport. 
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FIG. 1: Geometrical configuration of corrugated FS cylinder in the AF BZ (— 5 < k z < -5 with c = 2c), vortex direction (1) and 
quasiparticle momentum k. Here x, y = y, z is the cartesian crystal coordinate system with a,b and c (n)-axes. The cartesian 
coordinate system x,y,z is rotated by the angle 9 around y axis such that z ligns up with the vortex direction 1. Here cf> is the 
angle between heat current jQ and the field rotation (xz-) plane. Experimentally the perpendicular configuration <f> — 90° is 
used. Note that for uniaxial crystals the vortex polar angle 6 differs from the field polar angle Oh, this misallignment is given by 
Eq. IB II and shown in the inset of Fig. [5] The quasiparticle momentum k is described by cylindrical coordinates ($ = ck z ,(p) 
and vf _L FS denotes the quasiparticle velocity. The superfluid velocity is denoted by v s and a v , p are the vortex coordinates. 



III. SUPERFLUID VELOCITY IN THE LONDON LIMIT 

The superfluid velocity field v s (r) is obtained from the field distribution h(r) according to Eq. (J2J. In the London 
limit (Xi £i; i=a,c with |j=a, _L=c) the latter is determined by the equation 

h + AjjV x (V x h) + (Ai - A|)V x [n(n • V x h)] = z$ 0< 5(p) (4) 

where p=(x,y) are the cartesian coordinates of the plane perpendicular to the vortex direction 1 || z (Fig. ^| and 
$o = hc/2e is the flux quantum. In general h has components both parallel and perpendicular to 1 |27j . Here we 
neglect the latter and in the following assume h = h(x,y)l. The superfluid velocity and field distribution in the 
xy-plane are not circular around the vortex in the case of a uniaxial symmetry. However, we may apply a scale 
transformation p > p' or (x,y) — >(x',y') given by 

x' = — — = p 1 cos a v (5) 

A ll 

y' = — = p' sin a v (6) 

with p'=[x' 2 +y' 2 ]5 denoting the distance from the vortex center, a v the azimuthal angle around the vortex and 
9 = Z(l,n), see Fig.^ This leads to a transformed h(x',y') determined by 

h(x',y')- (d 2 x , +d*,)h(x',y') = ^-S(x',y') (7) 

Therefore in the scaled x',y' coordinate system one has again a circular vortex and field distribution. The solution of 
Eq. (7J is given by (K = Hankel function) 

h ^-4^ W) "4^ [1 v +0 - 12] (8) 

The scaled superfluid velocity field v' x =v x /Am, v'y—Vy/Xe for distances p' S> £/A is then obtained from Eq. J5| as 
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where A = (AiiAe) 2 , £ = (CllCe) 2 an d m — are defined in App.[5] The scaled p' coordinate is dimensionless. 

We may reintroduce length dimension but keep the circular vortex shape by defining 

p = \p' , v s =Av' s , (10) 

In the rotated coordinate system (x,y,z) of Fig. ^ with the z-axis alligned with the vortex direction 1, the scaled 
quasiparticle velocities defined in Add. 1X1 are transformed to 

' A^ 1 [w a (?9) cos 9 cos if + v c (d) sin 9] \ 
v f= _A e - 1 [u (^)sin^] (11) 

\ [— Va(&) sin cos (,5 + u c ($) cos 9] J 

Here «j(t?) = v^ii^) where t? is the polar angle for an ellipsoidal FS and ■& — ck z for the cylindrical FS of App. 
Furthermore (<£>,$) defines the direction of quasiparticle momentum on the FS and 9 is the vortex tilt angle in the 
field rotation plane which for convenience is chosen as the iz-plane. 

The dimensionless Doppler shift energy x (not to be confused with the cartesian coordinate) for a given vortex 
direction is then given by 

x = PJ? A Vs = ^rav F (jj, <p- 9) ■ v s {p, a v ) (12) 

with vf = Xv' F and v s defined by Eqs. (ff?Tfl . ffTTTfi and A denoting the SC gap amplitude. Using Eqs. (f^.lfTTJl we finally 
obtain for x = x(0, if; p, a v ) the explicit result 



- — ( (— -)v a (-&) sin if cos a v — (-r) [v a W cos 9 cos if + a^vJfl) sin 9] sina^l (13) 
4pL v A A J 



where y/a — v®/v® is the anisotropy of the Fermi velocity, £° = (Hv^/ttA) is the in-plane coherence length and Va,c(^) 
are given in App.^ Thus the Doppler shift energy of a quasiparticle depends on three sets of variables: The vortex 
direction (9) (or field direction 9h), see App.[BJ), the quasiparticle momentum coordinates ($, if) on the FS and the 
quasiparticle position in real space (p, a v ) with respect to the vortex center. The expression in Ea. i|13|) is the central 
quantity which determines the thermodynamic and transport properties in the vortex phase. In the special case for 
H || c (0 = 0) the DS simplifies to 

x= (0 )sm{<f-a v ) (14) 

4 p 

Because x now depends only on the angle difference if — a v , averaging over the DS for H || c involves one integration 
less than for general field direction. 



IV. AVERAGE OVER THE VORTEX COORDINATES 



The calculation of thermodynamic and transport coefficients and the quasiparticle DOS involves averaging over 
both vortex coordinates (real space) and quasiparticle velocities (momentum space) . It is instructive to calculate the 
vortex averaged DS energy (|x|)^ of quasiparticles with given k, i.e. as function field direction 9. Since the 

important contributions to 7(T,H) and «ii(T,H) come from the nodal regions with |A(k)|/A < \x\, the field-angle 
dependence of for k G nodal region gives already qualitative information on the behaviour of 7 and Ku as shown 

in this section. 

In the scaled coordinates x',y' spanning the plane perpendicular to the vortex direction 1(9) the field distribution 
h(x',y') and velocity distribution v s (x',y') are circularly symmetric for a single vortex. In the limit H H C 2 the 
average over the vortex lattice is then approximately replaced by the integration over a single rotationally symmetric 
vortex according to 

(A(p, a v )) v = — T =— / da v dppA(p,a v ) (15) 

TT{d z - dj.) Jo Jd c 
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FIG. 2: Vortex-averaged dimensionless Doppler shift (H = IT) as function of vortex direction 9 for various quasiparticle 
momenta k which are given in cylindrical coordinates as (i? = ck z ,0). Full line: (-§,0); dashed line: these k-vectors 

lie in the k z =7r/2c symmetry plane which is also the nodal AF zone-boundary plane for Ai 9 . dash-dotted line: (-| ,§); this 
k- vector lies in a non-symmetry plane which is the nodal plane for the last gap function in Table |l| The small splitting at 
6 = 90° is caused by the FS corrugation. The inset shows the misallignment between vortex (9) and field (Oh) directions for a 
mass anisotropy a=0.69 appropriate for UPd2Al3. In all following angular plots we use 9 as variable and neglect the difference 
between 9 and 9h- 



Here d c is a lower cutoff of the size of the coherence length £ = (^y^e) 2 and d is of the order of the inter- vortex 
distance. It is determined by assuming a square vortex lattice and replacing its (square) unit cell of area Fq by a 
circle of equal area with radius d requiring ird 2 = Fq=$q/H. This leads to 



with $o=(hc/2e) = 2.07T0 11 Tcm 2 and H = IT the magnetic length scale involved is do (IT) = 257 A. This means 
do = 3 when we approximate the BCS or Pippard coherence length £° by the a,c- averaged value £o ~ 85A given 
in Ref. ^| for UPd2Al3. This estimate has some uncertainty depending on the inclusion of Pauli limiting effects. 
The ratio <io/£a determines directly the size of the DS via the prefactor in Eq. (|14|l when p is expressed in units of do. 
Decreasing the ratio do/£fi at fixed H(do) increases the DS and hence the oscillation amplitudes. It is useful to check 
the consistency of this value in an independent way. Using the expression for H c2 in Eq. I|B2I) we may also write the 
magnetic length scale as 



do = V2C(^Y (17) 



With the experimental H° 2 = 3.2T [3 we obtain do (IT) = 215 A which is consistent with the previous value. 
The discrepancy may be due to Pauli limiting effects which reduce H C 2 from its purely orbital value. These 
estimates also give an insight in the validity range of the single vortex and DS approximation. For fields H <C H c2 
(but large enough to fulfill |a;|A T) we have £o ~ d c <C d, consequently we may approximate d c ~ in this 
limit. This means we commit a small error by extending the DS approximation to the core region where it is not valid. 

We now perform the average in Eq. 1)15(1 for A(p, a v ) — \x(9; $, (p\ p, a v )\, i.e. the absolute value for the DS energy 
in Eq. (|13J) for fixed momentum k (i?,y>) direction. In Fig. |3we compare (|jc|)„ for various quasiparticle momenta 
k. Keeping in mind that according to Eq. i|12|) \x\ vanishes when vp 1 v s and becomes maximal for || v s , the 
9 (or Oh) variation of {\x\y% in Fig. [21 can be qualitatively understood. It is the 8 variation of (\x\)^ for k £ nodal 
region which survives in the ZEDOS and transport coefficients from which conclusions on the positions of nodes of 
A(k) may be drawn. 



V. THERMODYNAMICS AND TRANSPORT IN THE VORTEX PHASE IN DS APPROXIMATION 



The low temperature transport and thermodynamics in unconventional SC are determined by the combined effect 
of impurity scattering and DS due to the supercurrents. Both effects may lead to a low energy residual DOS which 
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FIG. 3: Quasiparticle DOS for the Ai g gap function for H = (full line) and IT (dash-dotted line) and field || c [Oh = = 0). 
The DOS singularity is smeared out by the k-dependent DS. Inset shows the ZEDOS induced by the Doppler shift on an 
enlarged scale. 



influences specific heat and thermal conductivity. In the zero-field limit the theory is well developed 4\ . Our intention 
here is to study typical signatures of the nodes of A(k) in magnetothermal properties to draw conclusions on the gap 
structure. For this purpose the 'superclean limit' where the DS energy dominates the effect of scattering (|x|A ^> T) 
is the relevant one. 

Transport and thermodynamics of unconventinal SC in DS approximation has been developed by many authors over 
the years, following the pioneering work of Volovik [l^, we mention only a few of them here |30l l3lL l32l l33l 133 . l3fjl l3(| . 
For our purpose it is sufficient to have a summary of these results in the superclean limit in a concise form useful for 
numerical computation of the field-angle dependence of k^(H, T) and 7(H, T). In the zero-field case the quasiparticle 
energies of an unconventional superconductor are given by E k = [e k + A k ] s where A k is the nontrivial singlet gap 
function or A k = |dk| 2 with dk denoting the (unitary) triplet gap functions. The k (p) dependence of the gap 
function may be characterised by the form factor $(k)=A(k)/A where A(T) is the gap amplitude obtained from the 
solution of the gap equation. 

In the vortex phase the DS leads to an additional position (p,a v ) and field-angle (9) dependence of quasiparticle 
energies according to Eq. JTSJ. Defining E' k = E k /A we obtain for k (tf, ip): 

E k (9;p,a v ) = E k -x k (9;p,a v ) (18) 

Calculations of 7(6?) and Ku(9) therefore involves, in addition to the FS averaging present already in the zero-field 
case, the averaging over vortex coordinates as prescribed in Eq. (|15|) . 



A. Quasiparticle DOS and specific heat 

In the zero-field case (T — » 0) the quasiparticle DOS is given by 



N(E)/N = gi(E) with 9l (E) = Re g(E) 

g(E) = 9l (E)+ig 2 (E)= [ dS k — ^— T (19) 



■ [E 2 - A 



where dSk = j^dQ k = j- sin fld-dd<p or dSk =^pid(pdk z for the ellipsoidal and cylindrical FS case respectively. In the 
vortex phase one has to replace E by the Doppler shifted quasiparticle energies and form averages Gi^iE) = {gifl{E)) v 
over the vortex coordinates according the prescription of Eq. l|15fl . Then Gi^(E) and the field-angle dependent 
quasiparticle DOS are given explicitly as 



N{E 1 H 1 6)/N = G l {E) = (J dS u 



j= 2 \E k \(\El-Al\ + El-Al)i 
\K - A k | 



■4. |£ k |(|^-A k |-^ + A k )^ 
G, E] = ( / dS^ |^_ k A|| ), ^ 
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FIG. 4: Effective quasiparticle lifetime (Eq. 1241 as function of energy in Born approximation (B) and unitary (U) case for 
fields parallel {0—0) and perpendicular (#=90°) to c-axis (H = IT). The lifetime and its angular variation are much larger in 
the Born case. As shown in the inset (So = scattering phase shift) this leads to a much larger thermal conductivity oscillation 



amplitude 



; (90°)]t=o of up to 27% in the Born(B) case, whereas it is only ~ 1.2% in the unitary(U) case. The 



amplitude changes sign at S — 0.83(1-) and is maximal at <5o 



0.91|. 



where E^ is given by Eq. I|18|) and the vortex average (..)„ is defined in Eq. I|15|l . Note that N(E,H,0) also depends on 
the field strength H via the DS, this variable is sometimes not written explicitly. From Eq. (|2U[1 the field-angle and 
temperature dependence of the specific heat C(T,H,#) may be obtained in the usual way. Defining the linear specific 
heat coefficient as 7(T) = C/T and using 7„ = ^f-^o we obtain after a variable substitution (e = E/T) 



C(T,H,9) 7(T,J?,< 



In 



de- 



2n 2 N J cosh 2 (e/2) 



N(Te, H, I 



(21) 



Together Eqs. (|18|l . (|20|l and (|21|l allow to calculate the field-angle dependence of the specific heat in the vortex phase. 
These equations are valid within the DS approximation for the superclean limit for any anisotropic gap function A(k). 

Altogether a fivefold integration over momenta (i?, (p), vortex coordinates (p, a v ) and energy e has to be performed 
in general. For T — > one needs only the ZEDOS N(0,H,6>) and only four integrations are left. At this stage we have 
to proceed with numerical calculations to give definite quantitative predictions. Approximate analytical evaluations 
usually give only angle dependences but not the absolute magnitude of the DS effect. 



B. Magnetothermal conductivity 



It is an important advantage of the DS method as compared to the more advanced semiclassical methods |3jj, |38| 
that it provides expressions for both specific heat and thermal conductivity, whereas the latter method sofar can 
only be used for the DOS and specific heat. In the zero-field case the thermal conductivity in the SC phase can be 
calculated within the linear response approach of Ambegaokar and Griffin |39| . Using the DS approximation this has 
later been extended to the vortex phase [3(], Hjjl ^3 . 

In the normal state one has for the thermal conductivity 



2n 2 



N T ((Vi 



IFS 



T 



(22) 



and for the anisotropic FS case the 



where To = is the quasiparticle life time. In the isotropic case 
ratios ((v z ) 2 } fs/ ) 2 }fs are given in App.lAlbv Eq. (|A8|) . 

In the presence of a SC gap the new energy scale introduced by A leads to an energy dependent effective life time 
of quasiparticles in the vortex state which is given by Q 



r(E) 

TO 



= r(E) 

X(E) 
Y(E) 



X(E) 2 +Y(Ef 



Gi(E)X(E) + G 2 (E)Y(E) 
cos 2 6 + sin 2 S [Gl(E) - G 2 2 (E)] 

2 



(23) 



2 sin 2 6 G 1 (E)G 2 (E) 
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where Gi(E) and 62(E) have been defined in Eq. (|20|l . Furthermore So is the (isotropic) scattering phase shift which 
lies in the intervall [O,^]. In the limiting cases one has: 

So = : f(E) = (Born limit) 

8o = l : r{E) = Gl{E) G ^ [E) (unitary limit) (24) 

The low energy behaviour of f (E) in the Born and unitary limit is quite different as seen in Fig. therefore the low 
field behaviour of the thermal conductivity is very sensitive to the size of the scattering phase shift as discussed later 
(Fig.©. 

Thermal conductivity in the SC state for zero field involves a quasiparticle energy integration and FS momentum 
averaging This is a special case of the magnetothermal conductivity in DS approximation which we discuss here. 
The latter is obtained in the same spirit as the field-angle dependent j(T, H)-value: The SC quasiparticle energies are 
replaced by their Doppler shifted values according to Eq. i|18H and an additional averaging over the vortex coordinates 
has to be performed. Then we obtain the final result (e = E/T) 

«»(T,H) 3 f°° dee 2 . ((v 2 k K(Te, k, r))) Fsy 



«S(T) ^ 2 J cosh 2 (e/2) v ; (v&fs 

K(E,k,r) = -^-[E 2 - A(k) 2 ]?0(E 2 - A(k) 2 ) (25) 
\E\ 

where is the Heaviside function and E = AEk(6, p,a v ) whith E^ again given by Eq. (JTSJ. Note that i=x,y, z 
refer to the fixed crystal coordinate system, although for brevity we will later use the conventional notation n yy etc. 
without the tilde. The above expression for Ku(T) is on the same level of approximation as Eq. I|21|) for "f(T, H). 
It involves a five-fold integral due the FS averaging ($,¥>), vortex (v) averaging [p, a v ) and energy (e) integration. 
Finally we note that the expression for the effective lifetime in Eq. (|24|l is perturbative with respect to T. Since we 
consider only the superclean limit where A|x| ^S> T this is justified. 



Due to the cylindrical symmetry of both FS (Fig. ^ and gap functions the calculated thermal conductivity depends 
only on the relative angle 4> between field rotational (xz) plane and heat current jq. The dependence is caused 
by the factor u 2 k in the double average of Eq. [23] For jq parallel (k xx ) or perpendicular (n yy ) to the field rotation 
plane one has v 2 k = v a {k z ) 2 cos 2 {(p) and w 2 k = v a (k z ) 2 sin 2 (<^) respectively. For general angle one has to use 
= v a {k z ) 2 cos 2 (<^3 — 4>). Experimentally the perpendicular configuration with <j) — 90° has been used 14], therefore 
we focus on K yy (9,H). 



VI. APPLICATION TO UPd 2 Al 3 



The DS-based calculation scheme for magnetothermal properties described in detail before will now be applied to 
discuss recent field-angle resolved thermal conductivity measurements in UPd2Al3 01 . In this work it was established 
that the gap function of UPd2Al3 possesses a line node in the basal plane by a qualitative discussion of the experimental 
results. It was also argued that the experiment cannot distinguish between several possible gap functions proposed in 
PH which have different positions of the node line along the hexagonal axis. This was attributed to the UPd2Al3 FS 
geometry which is characterised by a dominating corrugated cylinder sheet oriented along c* . Later a further proposal 
implying a gap function with a node line lying in a non-symmetry plane was made [41| . Before discussing the results 
for these models of A(k) we summarise their basic symmetry properties and microscopic background. 



A. The symmetry properties of gap function candidates 

UPd2Al3 is the only HF superconductor whose microscopic mechanism of Cooper pair formation is known with 
some certainty. As mentioned in the introduction this is due to the partly itinerant and partly localised nature 
of 5f-electrons. The latter lead to a well defined magnetic exciton band seen in INS experiments. Most crucially 
the magnetic exciton at the AF wave vector Q is also seen in a strong coupling signature of the tunneling DOS 
of conduction electrons at an energy (~ 1 meV) that is identical to the INS results. This is strong evidence that 
magnetic excitons originating in CEF excitations of the local 5f subsystem mediate the Cooper pairing of itinerant 
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FIG. 5: The Ai s ZEDOS or specific heat 7(H,T=0) coefficient as function of field strength for 9 — 90° is shown as full curve. 
Typical \fH of the low field ZEDOS is observed. Thermal conductivity Kj^HjT^O) for 8 = 90° normalised to normal state 
value Ky y is shown for various scattering phase shifts So = 8/(ji/2) close to the unitary limit 8a — 1 (broken curves). Deviations 
from So = 1 lead to an increasingly prominent nonmonotonic behaviour of k yv (H) = K yy (H)/ n yy . The nonmonotonic field 
dependence was observed in [l4j|. For even smaller So approaching the Born case the initial drop is preserved but then h yy (H) 
becomes flat. 

TABLE I: Spin and orbital structure of the possible gap functions which are solutions of the Eliashberg equations for the 
dual model of UPd2Al3. The form factors of the anisotropic gap function A(fc z ) = A$(£; z ) are assumed to have cylindrical 
symmetry. The state in the last row |41| has a hybrid gap function since $(A: Z ) = 2 cos 2 ck z -l. The line nodes k z — ±7r/4c are 
due to fine tuning of the amplitudes of two fully summetric basis functions. 
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5f electrons. This mechanism has been investigated both in weak coupling 20j and in a strong-coupling Eliashberg 
approach |2l| . In the latter a model for the effective interaction based on magnetic exciton exchange with Ising 
type spin space symmetry was proposed. This breaks rotational symmetry in the (pseudo-) spin space of conduction 
electrons in a maximal way, therefore the pair states \x) have to be classified according to equal spin pairing (ESP) 
and opposite spin pairing (OSP) states characterised by a spin projection factor p =(|j \<J z cr z \ H) rather than in 
terms of singlet and triplet pairs. It was found [2l| that three of these states (type I-III) have a finite T c . The largest 
T c belongs to two degenerate OSP states of opposite parity. These states together with their orbital dependence and 
symmetry classification are tabulated in Table [I] In addition we have included a hybrid gap function (last row) of 
type IV proposed in |4l| consisting of a superposition of two inequivalent fully symmetric T>Qh representations. Due 
to this fact its nodal lines are lying in non-symmetry planes. This is rather uncommon feature and not observed in 
any unconventional SC sofar. In addition this gap function does not appear as a possible solution of the Eliashberg 
equations in the model of j^. Nevertheless we include it in the present discussion because it was proposed as a 
candidate in [4l| . 

B. Results of numerical calculations 

In the following we discuss the numerical results using the above analysis for DOS, specific heat and thermal 
conductivity. We will use the model parameters 7 C = 0.8 for the corrugated cylindrical FS and a = 0.69 for the 
anisotropy of the Fermi velocity appropriate for UPd2Al3. These parameters are taken as independent but in 
principle they are related via Eq. I|A6() . The impurity scattering phase shift So which is the only free or unknown 
parameter in the theory is mostly taken to be close to its unitary limit 7r/2 which is commonly assumed for 
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FIG. 6: Dependence of ZEDOS or specific heat 7(H)- coefficient on vortex direction 9 at different field strengths. Qualitative 
behaviour for all type I-III gap functions is equivalent, this also holds for K yy (9,H). 



HF compounds. We choose a value So = 0.9(7r/2) for which the oscillation amplitude of K yy (9,H) is close to 
its maximum. However we also study the thermal conductivity for more general <5o- Because we consider the 
superclean limit the ZEDOS is not influenced by the choice of Sq. We will discuss the typical results for the various 
gap functions in Table [I] but will not present an exhaustive overview of the results. The intention is rather to 
investigate whether the classification of gap functions introduced in 0] is justified. The four gap function examples 
in Table [i] correspond to the different types I- IV of nodal structure whose field-angle dependence of specific heat 
and thermal conductivity has been qualitatively discussed already in [T^j. Using the theory outlined in the previ- 
ous sections we can now perform detailed numerical calculations and check the conjectures given in |l4j quantitatively. 

The important effect of the Doppler shift of quasiparticle energies shown in Fig. [21 is the appearance of a 
non- vanishing ZEDOS (Fig. [3J) . Since the DS depends strongly on the quasiparticle momentum k for a given field 
direction the ZEDOS, which is dominated by quasiparticles in the nodal regions, will exhibit pronounced field-angle 
dependence in addition to its dependence on field strength. To simplify the discussions in the following we do not 
distinguish any more between field (Oh) and vortex (0) directions since they are very close for the present value of 
a (see inset of Fig. [2J|. The field dependence of ZEDOS or specific heat 7-coefficicnt is shown by the full curve in 
Fig. for the Ai g gap function at 8 — 90°. It exhibits the typical ~ \fii- behaviour for nodal gap functions which 
is due to the DS shifted continuum states in the inter- vortex region as first predicted by Volovik [13j. This is in 
contrast to the ~ H behaviour of the specific heat coefficient in isotropic superconductors for H <C H c2 which is due 
to the quasi-bound states in the vortex cores. The field dependence of K yy for 0=90° is shown in the same Fig. [S] for 
various scattering phase shifts. The behaviour is closer to linear H-dcpcndcncc. If the scattering phase shift deviates 
from the unitary limit just a few percent, then immediately a nonmonotonic low-field behaviour with an associated 
minimum in K yy (H) appears as is obvious from Fig. [5J It is caused by the low energy behaviour in the effective 
quasiparticle life time shown in Fig.^J Such nonmonotonic behaviour has indeed been found experimentally in [T^ |. 

The field-angle dependence of the ZEDOS or specific heat 7- value as shown in Figgis rather directly determined 
by that of the DS of nodal quasiparticles as may be seen by a comparison with Fig. [21 keeping in mind the fact that 
it becomes large for a field angle when quasiparticles with all k - vectors are Doppler shifted as is the case for 6=0. 
The maximum of the oscillation observed for this angle increases with \ H behaviour as mentioned before. For fields 
reasonably well below H C 2 the oscillation amplitude is of the order of several per cent of the normal state DOS No 
or normal state 7™ (Fig.[SJ|. In this figure we also show a comparison of angular dependence for the first three order 
parameters (type I-III) in Table [i] which all have node lines parallel to the hexagonal plane but at different values 
of k 2 . In Ref. hj it was suggested with qualitative arguments that the angle dependences should be similar in the 
three cases. The reason is that the Doppler shift is determined by the product v s • v> and vp is always parallel to 
the hexagonal plane for the possible node line positions k z = ±f ,0 of type I-III. Since the latter are reflection planes 
the Fermi velocity perpendicular to these planes has to vanish as is obvious from the corrugated cylinder FS of Fig.^ 
This qualitative expectation is indeed confirmed by the numerical results for type I-III gap functions in Fig. [H] which 
shows the same type of 9 dependence. The small difference of the amplitude is due to the different size of the Fermi 
velocities at the above node line positions (Eq. IA6J) due to FS corrugation. The present fully numerical treatment 
of the DS theory allows also to calculate the temperature dependence of the oscillation amplitude of 7(0, T); it is 
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FIG. 7: Temperature dependence of the normalised field-angular oscillation amplitude for specific heat or ZEDOS at different 
fields. The amplitude is defined by Sj — 8N — [N(0) ~ N(90°)]t/[N(0) - 2V(90°)] . 



shown in Fig. [7| The amplitude first decreases with T 2 behaviour and then decreases rapidly. For UPd2Al 3 one has 
2A a „/kT c = 5.6 0, therefore (T/T c ) = 2.8(T/A aw ). Here A av should be interpreted as a FS-sheet and momentum 
averaged gap value. If we only consider the cylindrical FS sheet, then for all of our gap models (I-III) we have 
A av =A/\/2 which means (T/A)=0.25(T/T C ). Comparison with Fig. [7| shows that the T=0 oscillation amplitude 
has dropped to 10-15% of its original value when the temperature has increased to 20% of T c . This illustrates the 
neccessity of having T<cT c if one wants to observe DS induced angular oscillations of magnetothermal quantities. 

Similar observations can be made for the thermal conductivity na{9) for type I-III gap functions as shown in Fig. [HI 
where we discuss K yy (9). This component corresponds to the experimental configuration with heat current perpen- 
dicular to the plane of field rotation. The calculation has been done for the nearly unitary phase shift 8q=0.9(tt/2) 
where the oscillation amplitude is close to its (positive) maximum as shown in the inset of Fig. 0] For this 6 value 
K yy (90° , H ) has a minimum around H m ~ 0.4T (Fig. |3J). For fields above H m again the maximum appears at 9=0 
when the field is perpendicular to the nodal plane. However for fields below H m a sign change of the amplitude takes 
place and therefore the maximum appears for field direction lying within the nodal plane. This is due to the increased 
influence of life time angular dependence at small fields. The shape of angular oscillations, their field sequence (Fig. [5] 
left panel) and the minimum of n yy (90°, H) in Fig. [5] correspond well to the experimental observations in [l^ |. 

The relative differences in the oscillation amplitude for type I-III gap functions caused by FS corrugation are 
somewhat larger as in the case of ZEDOS, however their qualitative behaviour is again indistinguishable and is not 
shown here. The absolute oscillation amplitude of the order of per cent of the normal state conductivity at IT is 
smaller than for ZEDOS (Fig. |BJ). In the present superclean limit the latter is independent of the scattering phase 
shift, i.e. unitary or Born limit. In contrast the absolute oscillation amplitute of k(9,H) depends strongly on the 
phase shift 5o via the pronounced So dependence of the effective life time as seen in Fig. 0] As the inset of Fig. 0] 
shows, the oscillation amplitude for n yy (9,H) is very large (with opposite sign) in the Born limit and comparatively 
small (of order per cent) in the (nearly) unitary limit which we have assumed here for UPd2Ai3. A similar calculation 
for n xx {9, H), i.e. heat current parallel to the field rotation plane shows that the oscillation amplitude is always 
positve and again much larger for Born as compared to nearly unitary scattering. The calculated absolute magnitude 
of the T=0 thermal conductivity n yy (9, H) is smaller by a factor of three compared to the experimental value at the 
lowest temperature measured. In addition a similar T-dependendence as for ZEDOS (Fig. [TJ) would lead to a further 
strong reduction. These discrepancies may be linked with our insufficient model for the energy and field dependence 
of the quasiparticlc life time. They cannot presently be resolved without experimental results on the T-dependence 
of the oscillation amplitudes. For even larger larger fields (H > 2.5T) than in Fig. |S|an additional sign change in the 
oscillation was observed. In this regime the DS approximation breaks down due to vortex overlap and the sign change 
was rather attributed to the influence of the uniaxial H C 2-anisotropy [l4j . 

Since angular oscillations are not observed when the field is rotated in the hexagonal plane |14| (and heat current 
parallel to c-axis) one may conclude that the node line of the gap is indeed parallel to the hexagonal plane as for all 
type I-IV gap functions discussed here. Our calculations prove quantitatively that the observed K TO (0)-dependence 
cannot distinguish between the type I-III cases since the 9 dependence for these gap functions is very similar as has 
already been suggested in 0] on qualitative grounds. The reason has been discussed above in the context of ZEDOS 
(Fig. [SJ oscillations. 
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FIG. 8: Left panel: Thermal conductivity as function of field angle for type II gap function at various field strengths. This 
behaviour, notably the sign change of the oscillation amplitude with increasing field corresponds to experimental observations 
for fields H < 2T in fhH ]. Type I, III gap functions (not shown here) exhibit very similar oscillations as shown in the case 
of ZEDOS (Fig. Right panel: Thermal conductivity for the hybrid (type IV) gap function A(k)= Acos2cfc z . A sharp 
minimum at intermediate 9 ~ 51° due to off-symmetry nodal plane appears. This behaviour is in contradiction to experimental 
observations [T^l. 



The situation however is different for the hybrid Ai g <g>A' lg gap function (type IV) which has node lines in non- 
symmetry planes. The corresponding angle dependence of the thermal conductivity k xx (8) is shown in the right panel 
of Fig. |SJ It has a completely different appearance from those for the type I-III gap functions. Firstly, they do not 
drop to a small value for field angle 9 = ±90°, instead they are even larger than for 9—0. Secondly a pronounced and 
sharp minimum appears at an intermediate angle 9 ~ 51°. This is due to the existence of a non-symmetry nodal plane 
for the Ai g £g)A'ig order parameter. In such case both v x , y (^z) and v z (k z ) are nonzero and they give contributions to 
the DS in Eq. 112fl which decrease or increase as function of 9 respectively, leading to the minimum at an intermediate 
value that survives in the averaged quantities like thermal conductivity. The minimum cannot be discussed away by 
including the contribution from other FS parts. Firstly it is known from LDA calculations 0] that the cylinder FS 
gives one of the dominant contributions to the total DOS, secondly the other, e.g. ellipisoidal sheets also have a finite 
v 2 at a non-symmetry nodal plane and therefore would give a similar behaviour. Consequently, if the gap function is 
of the type IV the kink- like minimum at intermediate 9 has to be present in k xx {9). However, as mentioned above, 
the experiments Q show that it behaves very much as expected for the type I-III gap functions discussed before 
(Fig. [SI left panel). No trace of a minimum or only depression at intermediate angles has been found. Therefore one 
has to conclude that the experiments in 0] rule out the type IV gap function A(k)=A cos(2ck z ) proposed in for 
UPd2Al3. It was already suggested before that in any case this is an unlikely candidate because it is a hybrid gap 
function with nodes at non-symmetry planes, which has never been found in any other unconventional superconductor. 



VII. SUMMARY AND OUTLOOK 



In this work we have investigated the ZEDOS, specific heat and magnetothermal transport properties of supercon- 
ducting UPd2Al3 in the vortex phase using the Doppler shift approximation. Our intention was to give a quantitative 
basis to the qualitative discussion of possible gap function symmetries presented together with the experimental 
results of [l4j . We have given a coherent presentation of the known Doppler shift analysis and expressions and 
performed the evaluation of physical quantities in a fully numerical approach which allows quantitative predictions. 
The DS approach is oversimplified in the sense that it does not correctly account for the vortex core contributions 
and the effect of vortex overlap on approaching H C 2. Therefore one is limited to fields reasonably well below H C 2 
but still large enough for the superclean limit to hold. On the other hand it has the great advantage as compared 
to more advanced semiclassical methods that transport quantities and not only the quasiparticle DOS are easily 
accessible within a linear response treatment. We have made predictions for the candidate type I-IV gap functions of 
UPd2Al3 which have all line nodes parallel to the hexagonal plane but with different multiplicity and position along 
k 2 . Our quantitative analysis fully confirms the qualitative conjectures drawn in 14]. Only gap functions which 
have nodal lines in symmetry planes containing the BZ boundary or center are compatible with experimental results, 
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those with nodal lines in off- symmetry planes can clearly be ruled out. It was suggested that the nonmonotonic 
behaviour of the low-field thermal conductivity may be caused by the sensitivity of the effective life time to devi- 
ations from unitary scattering. This is also connected with the sign change of the oscillation amplitude for small fields. 

The angular- resolved thermal conductivity thus proves that a node line of A(k) must be present in a hexagonal 
symmetry plane, however by itself it cannot distinguish between the possible gap functions of type I-III. As explained 
in [l4j one needs additional information: Inelastic neutron scattering |17| requires the translation symmetry A(k + Q) 
= -A(k) which would exclude A'i„, also this gap function is disfavored in Eliashberg theory [211 . Furthermore naive 
interpretation of Knight shift results |42| advocates for an even parity gap function, although a comparitive analysis 
for the Aip and Ai u states has not been performed [2l| yet. This finally led to the suggestion 0] that the Ai ff gap 
function in the second row of Table ID is the proper one for UPd2Al 3 . 

Presently we have restricted ourselves to gap functions A(i9) that have cylindrical symmetry. Our numerical method 
can straightforwardly be applied to fully anisotropic gap functions A(t?,y>). The calculation then provides us with 
a direct mapping of A($, (p) — > j(0h,4>h) or ^u^Ht^h) between the momentum dependence of gap functions and 
the field-angle dependence of physical quantities in their respective 2D domains. Of course this mapping has no 
unique inverse, but nevertheless comparison with the experimental 9jj, 4>h- dependence may provide important clues 
on the anisotropy character of the gap functions. In addition, an extension of the present quantitative theory to 
different types of FS sheets like FS ellipsoids, corrugated tight binding with FS nesting features etc. is easily possible. 
Finally we note that the numerical approach to the DS theory may in principle be generalised beyond the perturbative 
treatment of scattering, i.e. from the superclean to the clean limit with A|x| < T <C A, although this would likely 
mean a much larger computational effort. 



APPENDIX A: QUASIPARTICLE VELOCITIES FOR CORRUGATED CYLINDRICAL FS 

In this appendix we define the geometric features and quasiparticle properties of the corrugated cylindrical Fermi 
surface (FS) which is necessary for the calculation of the Doppler shift energies. Specifically we give the quasipar- 
ticle velocities in terms cylindrical coordinates. The corrugated cylinder is the most prominent heavy FS sheet in 
UPd2Al3 obtained in LDA calculations 0,0], dual model calculations ^|| and also from dHvA experiments (lif . 
The latter show it has also among the heaviest quasiparticle masses. In the AF BZ (—5 < k z < ? ) with c = 2c 
appropriate for UPd2Al3 it can be modeled as 

e(k) = e a (k)-2i c coscfc z , e a (k) = ^—(k 2 x + k 2 ) (Al) 

2m a y 

where the first part is the parabolic ab-plane dispersion and the second part is the tight-binding like dispersion which 
determines the FS corrugation along c. The diameter of the corrugated cylinder is given by 

k a F (k*) = k° F [l + — co S (ck z )]i ; fc°=(^£)* (A2) 

We introduce the FS corrugation parameter by 

lc = WT) ' ^ (A3) 

where 7 C is the ratio of FS cross sectional areas at the (AF) zone center and zone boundary respectively. It is given 
by 

The quasiparticle velocities = j-Vek in cylindrical coordinates are given by 

v x {k z ,0) = v°v a (k z ) cos ip 

V y (k z ,tp) = v"v a (k z ) simp 

v z {k z ) = v°v c (k z ) 

v° a = hk F /m a (A5) 
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v° c = 2-Kt c /hQ 

Va(k z ) = [1 + Acoscfc z ]3 
v c (k z ) = sincfc 2 

the anisotropy ratio of the Fermi velocity of quasiparticles is then given by (Q = ?): 

a= (4) 2 = ^(S) 2 A 2 or a =^(^!) 2 (2^I) 2 (A6) 

The anisotropy ratio a and corrugation factor j c are the parameters which determine the Doppler shift energy for 
the present FS. The relation between a and j c in Eq. (|A6I) is valid only for parabolic in-plane dispersion. It is better 
to assume these parameters as independent and take their values from experiment, keeping in mind that a — > for 
7c ^1. 

Finally the FS averages over quasiparticle velocities are given by 
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APPENDIX B: CRITICAL FIELDS, COHERENCE LENGTH PENETRATION DEPTH IN UNIAXIAL 

SUPERCONDUCTORS 

The theory of critical fields H c i and H C 2 in superconductors with uniaxial symmetry like and T)eh was given in 
[27L l2Sj (see also 0|) on the basis of Ginzburg-Landau theory for a single component SC order parameter. Here we 
give a summary of relations derived in these references which are important in our context of Doppler shift calculations. 
In uniaxial geometry shown in Fig. ^ coherence length and penetration depth are different for fields directed along a 
or c crystal axes. As a consequence, for intermediate polar field angle Oh (with respect to c) in the range < Oh < § 
the field (Oh) and vortex (0) directions are not the same. This is an essential difference to the isotropic case where 
field and vortices are alligned. For uniaxial effective masses m a , m c these angles are related by 



1 . „ x „ _i , . m a ,v° 



= tan _1 (-tan(V) ; H = tan _1 (a tan (9) ; a = — = (^f (Bl) 
a m c v y l 

where we assume the convention — ^ < 0,0h < ■§ for the vortex and field direction. The function 0(0 h) is plotted 
in the inset of Fig. [21 for the effective mass anisotropy a — 0.69 which is the appropriate average value for UPd2Al3. 
For cn-values only moderately different from one, and Oh are rather close. However, for a <C 1 the vortices (0) are 
'pinned' along the c-axis and they are lagging behind the field direction (Oh ) when it is continuously changed from c 
to a. The direction dependence of critical fields is given by 



H c2 (0 H ) = ^ ; $0 



2< 2 ' 2e 



. [(H^) 2 cos + (H c cl ) 2 sin 0] 
Where £ is given in Eq. IB5I and the uniaxial H°{ c are obtained from 



H ^-^^r ' — 27r with (B3) 

The uniaxial Ginzburg-Landau coherence lengths £ , c and penetration depths X a , c are given by (i = a,c) 

= and A 2 .-^- (B4) 
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with r=l-T/T c . The FS averages for the corrugated cylindrical FS are derived in Add. 1X1 The effective field-angle (Oh) 
dependent coherence length £ in Eq. (|_B2(I . penetration depth A and mass m are given by (An = A a , ^ = £ a , mil = m a ) 



e=te|6)», A=(Ai|Afl)» > m=(mj,m e )^ (B5) 



where we used the definitions 



Q cos 2 6> + sin 2 

8Trn s e 



2 

\2 _ \2„„„2/) i \2 • 2 /) _ m G C 



A 2 = A 2 cos 2 9 + A 2 sin 2 = £^ (B6) 



= m a cos 2 + m c sin 2 9 

(B7) 

Note that 9 is the vortex angle with respect to the c-axis which is related to the field angle 9h via Eq. i|Bl|) . From 
the above equations we derive the relations 

^- = (-!)' = a 3 [sin 2 (9 + a cos 2 6>]~ 3 
A Xe 

= (hly = a -i[ s [n 2 9 + aco S 2 9]i (B8) 
A A,, 

which enter directly the expressions for the Doppler shift energies in Eq. I|13(l . 
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